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ABSTRACT 

We study the internal response of a galaxy to an unbound encounter and present a survey of or- 
bital parameters covering typical encounters in different galactic environments. Overall, we conclude 
that relatively weak encounters by low-mass interloping galaxies can cause observable distortions in the 
primaries. The resulting asymmetries may persist long after the interloper is evident. 

We focus our attention on the production of structure in dark halos and in cluster ellipticals. Any 
distortion produced in a dark halo can distort the embedded stellar disk, possibly leading to the formation 
of lopsided and warped disks. We show that distant encounters with pericenters in the outer regions of 
a halo can excite strong and persistent features in the inner regions. Features excited in an elliptical 
are directly observable and we predict that asymmetries in the morphologies of these systems can be 
produced by relatively small perturbers. For example, a fly-by on an orbit with pericenter equal to the 
half-mass radius of the primary system and velocity of 200 km/s (a value typical for groups) can result 
in a significant dipole distortion for perturbers with mass as small as 5% of primary's mass. 

We use these detailed results to predict the distribution of the A parameter defined by Abraham et 
al. (sensitive to lopsidedness) and the shift between the center of mass of the primary system and the 
position of the peak of density for a range of environments. We find that high-density, low- velocity 
dispersion environments are more likely to host galaxies with significant asymmetries. Our distribution 
for the A parameter is in good agreement with the range spanned by the observed values for local galaxy 
clusters and for distant galaxies in the Medium Deep Survey and in the Hubble Deep Field. Assuming 
that primordial galaxies were located in dense environments with previrialized low velocity dispersions, 
our conclusions are consistent with the observational results showing a systematic trend for galaxies at 
larger redshifts to be more asymmetric than local galaxies. 

Finally, we propose a generalized asymmetry parameter A(r) which provides detailed information on 
the radial structure of the asymmetry produced by the mechanism explored in our work. 

Subject headings: celestial mechanics, stellar dynamics — galaxies: structure — galaxies:kinematics and 
dynamics — galaxies:clusters:general — galaxies:halo — galaxiesunteractions 



1. INTRODUCTION 

Galaxy interactions are likely to play a key role in de- 
termining the morphology and the structural properties 
of galaxies and in driving evolution at all epochs and in 
widely- varying environments. The importance of envi- 
ronmental conditions on the properties of populations of 
galaxies has been recognized for a long time (Hubble & 
Humason 1931) and observational studies have now firmly 
established relationships between the relative abundance 
of galaxies of different morphological types and the den- 
sity of the environment (see e.g. Dressier 1980, Dressier 
et al. 1997), their location in clusters (Oemler 1974) as 
well as the evolution of the morphology of galaxies in clus- 
ters at different epochs (Butcher & Oemler 1978, 1984). 
In particular, the ratio of early to late-type galaxies has 
been shown to be an increasing function of density and a 
decreasing function of the distance from the center of a 
cluster. In addition, HST observations show that the rela- 
tive number of spirals to elliptical and SO galaxies tend to 
increase in distant clusters, pointing to a time evolution of 
the morphology of galaxies from spirals to ellipticals (see 
e.g. Dressier et al. 1994a, 1994b, Couch et al. 1998). Al- 
together, these results suggest that environmental distur- 
bances due to interactions with other galaxies and with 
a cluster tidal field can lead to dramatic changes in the 



structure of galaxies. 

Simulations have explored strong encounters and merg- 
ers. Our recent work suggests that relatively weak en- 
counters can give rise to significant asymmetries as well. 
To study weaker galaxy interactions and quantify their 
importance, we compute the response of a spherical stel- 
lar system to the perturbation induced by another system 
during a fly-by, focusing our attention on the morphology 
and the persistence of the distortions produced during the 
encounter. The range of orbital parameters considered en- 
compasses the typical environmental conditions for field 
galaxies, for galaxies in groups and in clusters. Our per- 
turbative solution for these distortions provides quantita- 
tive and qualitative predictions in low to moderate am- 
plitude that are difficult to reach by simulation. This ap- 
proach nicely reveals the underlying dynamics without be- 
ing masked or swamped by n-body fluctuation noise. On 
the other hand, numerical feasibility restricts the complex- 
ity of possible geometries and the perturbative assumption 
limits its validity for large amplitudes. 

Our investigation is relevant for galactic dark halos or 
for spheroidal galaxies over a range of environments. We 
explore the possibility that the halo distortions produced 
during a fly-by could in turn lead to a significant and ob- 
servable distortion in an embedded stellar disk (see e.g. 
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Weinberg 1995, 1998b, Murali & Trcmaine 1998, Murali 
1998). In particular the deformations produced in a dark 
halo during an encounter could be one of the culprits for 
distorted morphologies like those observed in lopsided and 
warped galactic disks (sec e.g. Richtcr & Sancisi 1994, 
Zaritsky & Rix 1997, Rudnick & Rix 1998, Swatcrs ct al. 
1998, Haynes et al. 1998, Kornreich, Haynes & Lovelace 
1998). As recently shown by Weinberg (1995, 1998b), this 
could be the case for the Milky Way in which the deforma- 
tions induced in the dark halo by the Magellanic Clouds 
could be responsible for the observed warp in the Galac- 
tic disk. In a recent investigation, Reshetnikov & Combes 
(1998, 1999) have found that about 40% of a sample of 540 
galaxies have warped planes and their analysis revealed 
that warped galaxies tend to be located in denser environ- 
ments supporting a scenario in which interactions would 
play an important role in the formation of warps. 

Distortions of spheroidal galaxies are directly observ- 
able. Our results show that in low velocity dispersion 
and dense environments, where low-velocity close encoun- 
ters are more likely to occur, even low-mass perturbers 
can lead to significant asymmetries in the primary galaxy. 
Recent morphological studies of both local and very dis- 
tant galaxies have yielded comparisons of the frequencies 
of systems with distorted morphologies in different envi- 
ronments and at different distances (see e.g. Zepf & Whit- 
more 1993, Mendes de Oliveira & Hickson 1994, Abraham 
ct al. 1996a, 1996b, van den Bergh et al. 1996, Naim, Rat- 
natunga & Griffiths 1997, Brinchmann et al. 1998, Mar- 
leau & Simard 1998, Conselicc & Bershady 1998, Con- 
selice & Gallagher 1999). On the theoretical side, the 
efforts have been largely concentrated on numerical N- 
body simulations aimed at investigating the fate and the 
evolution of merging galaxies (see e.g. Barnes & Hern- 
quist 1992, Barnes 1998 and references therein) and the 
distortions and the morphological alterations produced in 
galaxies in clusters and groups by the ensemble of "tidal 
shocks" with other members (Moore et al. 1996, 1998). 
This mechanism can drive the evolution of spiral galaxies 
into spheroidal systems and it is a strong candidate expla- 
nation for the difference in the relative number of spiral to 
elliptical galaxies in clusters at different redshifts. 

The plan of the paper is the following. In §2 we out- 
line the method adopted for this investigation while the 
mathematical details are described in the Appendix. In 
§3 wc show the resiilts obtained for initial conditions typi- 
cal of galactic dark halos and their dependence on the or- 
bital parameters of the perturber. In §4 initial conditions 
relevant for spheroidal galaxies and a survey of orbital pa- 
rameters of the perturber typical of different environments 
are considered. Overall, our results quantify the relation- 
ship between the interaction rate for the environment and 
the probability of observing irregular morphologies, aid- 
ing interpretation of any trend in the fraction of asym- 
metric galaxies at different redshifts. This latter aspect 
is of particular interest as data on the morphology of dis- 
tant galaxies are now available from HST observations in 
the Hubble Deep Field and in the Medium Deep Survey 
projects (Abraham et al. 1996a, 1996b). In order to fa- 
cilitate the comparison of our theoretical results with ob- 
servational data, we will quantify the distortions produced 

1 Although this approach can treat a perturber with any density p 
by a point mass. 



by means of the asymmetry parameter A (Abraham et al. 
1996a) which has been determined for a number of local 
and distant galaxies. We will also introduce a generalized 
asymmetry parameter A{r) which will provide an impor- 
tant information on the radial structure of the asymmetry 
produced by the mechanism we have considered and which 
will allow to test the theory presented here. The main con- 
clusions are summarized and discussed in §5. 

2. METHOD 

This section outlines the perturbative solution of the 
coUisionless Boltzmann equation used to predict the dy- 
namical response of galaxies to interactions. Additional 
mathematical detail is presented in the Appendix. 

The linear perturbation theory used here has been ap- 
plied in the past both to investigate the coUisionless stabil- 
ity of stellar disks (Kalnajs 1977) and of spherical stellar 
systems (Polyachcnko & Shukhman 1981, Palmer & Pa- 
paloizou 1987, Bertin & Pegoraro 1989, Saha 1991, Wein- 
berg 1991, Bertin et al. 1994; see also Palmer 1994 and 
references therein) and to study the response of a galaxy to 
perturbations induced by satellite companions (Weinberg 
1989, 1995, 1998b). This method is numerically intensive 
but it avoids the problems of n-body simulations due to 
the finite number of particles which can introduce spuri- 
ous noise effects (see Weinberg 1998a) and it guarantees 
the resolution of resonances that lead to the excitation of 
patterns in the primary system. With too few particles, 
for example, noise can cause orbital energies and angular 
momenta to drift on a time scale shorter than the pattern 
speed of the global response of interest. This eliminates 
the possibility of observing this structure in the simulation. 
No doubt, real galax;ies are not smooth and inherent fluc- 
tuations may be important to their evolution. However, 
many millions of n-body particles are required to approach 
the estimates of natural amplitudes (Nelson & Tremaine 
1996, Weinberg 1998a). 

We compute the perturbation induced on a spherical 
stellar system by a point mass perturber ^ with mass m 
on a rectilinear trajectory with pericenter p. The response 
of the primary system will be determined by the simul- 
taneous solution of the linearized coUisionless Boltzmann 
and Poisson's equations, 

dh dhdHo _ dfodlh ^ , . 

m dw di di dw ' ^ ' 

V2$i =47rG/9i, (2) 

where the subscript denotes the equilibrium quantities, 
subscript 1 denotes the first order perturbation of the dis- 
tribution function /, the Hamiltonian H, and the potential 
The Boltzmann equation has been written in terms of 
action-angle variables I,w. The perturbed potential, $i, 
is the sum of the tidal potential due to the perturber, 
$p, and the gravitational potential of the response to the 
perturbation, ^respi and Pi = Pi + P\''^^ ■ The coupled 
Boltzmann-Poisson equations lead to an integro differen- 
tial equation for the potential. This can be solved by 
performing a Fourier transform in the angle variables, a 
Laplace transform in the time variable and expanding the 

le, there is no change in the global response if the interloper is replaced 



Perturbations of spherical stellar systems 



3 



perturbed quantities in terms of spherical harmonics and 
of a bi-orthogonal density-potential basis functions for the 
angular and the radial parts respectively. The integrodif- 
ferential equation is then reduced to an algebraic equation 
for the vector of coefficients of the expansion of the re- 
sponse potential and density, a, 



a = R(a + b). 



(3) 



The matrix R contains the time dependence of the per- 
turbation and the equilibrium properties of the perturbed 
system and b is the vector of the coefficients of the expan- 
sion of the external perturbing potential. 

It is clear from equation (||) that the self-consistent re- 
sponse of the perturbed system is determined both by the 
perturbation applied by the external perturber and by the 
reaction of the system to its own response to this pertur- 
bation. In general, an additional contribution to the to- 
tal response comes from excitation of the discrete damped 
modes of the primary system besides the effects of the per- 
turber and the self-gravity of the distorted primary. The 
differences in the mathematical structure of these two con- 
tributions are illustrated in the Appendix. The transient 
response due to the discrete damped modes can have a 
strong effect on the amplitude of the response (cf. §3) and 
it is particularly important when the modes are weakly 
damped with damping timescales longer than the fly-by 
timescale. In this case, the discrete damped modes pro- 
duce long-lived features which completely determine the 
amplitude and the structure of the response long after the 
pericenter passage of the perturber. Once the above alge- 
braic equation has been solved, the calculation of response 
potential and density is straightforward 
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where di and Ui are the bi-orthogonal density-potential 
basis function and Yim{0, (j)) are the spherical harmonics. 

3. PERTURBATION OF DARK HALOS 

Features excited in the halo can give rise to visible dis- 
tortions in an embedded disk even though dark halo per- 
turbations are not directly observable. While several other 
processes could be responsible for the observed distorted 
morphologies of stellar disks (see e.g. Sellwood & Merritt 
1994, Levine & Sparke 1998, Bertin & Mark 1980, Nelson 
& Tremaine 1995), the role played by a distorted galactic 
dark halo in exciting peculiar features in disks has been 
advanced for the case of an interaction of a galaxy with a 
satellite companion by Weinberg (1995,1998b) and appHed 

^the logarithm of the ratio of the total, Rt, to core radius, Rc 



to the case of the Milky Way-Magellanic Cloud system. 
Halo deformations leading to morphological and kinemat- 
ical lopsidedness has been invoked by Swaters et al. (1998) 
in an analysis of two kincmatically lopsided spiral galax- 
ies (DDO 9, NGC 4395). The ubiquity of MageUanic-class 
dwarf companions of normal spirals (e.g. Zaritsky & White 
1994 and Zaritsky et al. 1997) suggests that such inter- 
actions may be important sources of structure. Unbound 
encounters can have similar effect as we will illustrate here. 

Our standard dark halo model is an isotropic King 
model with a dimensionless central potential Wq = 3.0 
(concentration ^ c = 0.67) a total mass M ~ Q x IQ^^Mq 
and a total radius Rt = 200 kpc. These values are ap- 
propriate for the Galactic halo (see e.g. Kochanek 1996) 
and together with a standard exponential disk results in 
a fair fit to the observed rotation curve. Given the pri- 
mary model, the perturber's velocity V and pericenter p 
determine the structure of the response to the perturba- 
tion. The amplitude of the response scales with the mass 
of the perturber. 

Table 1 summarizes the values of V and p explored here 
along with the corresponding values of the ratio of the 
pericenter to the half-mass radius of the primary system, 
p/ Rh, and the ratio of the characteristic frequency of the 
motion of the perturber to the circular frequency at the 
edge of the primary system, Vl = V/p^J R^ /GM . All the 
cases considered here correspond to fly-bys in which part 
of the orbit passes through the halo; external encounters 
are likely to produce weaker distortions and will be consid- 
ered later in §3.3. The values of V considered are chosen 
to cover a typical range of relative encounter velocities in 
the field, groups and clusters. In order to clearly assess the 
importance of weakly damped modes, we consider the evo- 
lution of the response both with and without their effect 
for all the cases listed in Table 1. 

Table 1 

Orbital parameters for the fly-by 
encounters with a low-concentration 
primary system 



V (km/s) 


p (kpc) 


P/Rh 


VL 


200 


53.6 


1.0 


6.45 


500 


53.6 


1.0 


16.1 


1000 


53.6 


1.0 


32.3 


200 


107.2 


2.0 


3.23 


500 


107.2 


2.0 


8.1 


1000 


107.2 


2.0 


16.1 



For cases that include the weakly damped modes, we 
must first determine their natural frequencies. We used 
the method adopted in Weinberg (1994): the determinant 
of the dispersion matrix D(s) (see Appendix for its defi- 
nition) is evaluated on a grid of values of s in the upper 
{Re{s) > 0) complex plane where the determinant is de- 
fined and then rational functions are used to perform the 
analytical continuation to the lower complex plane where 
the zeros of the determinant corresponding to damped 
modes are located. Rational function fits of the elements 
of the dispersion matrix have been adopted also for the 
calculation of the residue of the inverse of the dispersion 
matrix necessary to include the effects of damped modes 
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in the response coefficients (see eq. Alt). The characteris- 
tic damping times for the most weakly damped modes, r/, 
are Ti/tdyn — 150 and T^jtdyn — 30 for the dipole and the 
quadrupole mode respectively, where t^yn is the dynam- 
ical time of the primary system defined at the half-mass 
radius, tdyn = tt^/B^J^GM. 

3.1. Importance of damped modes 

Figure 1 shows the evolution of the energy in the pertur- 
bation, Ep, as a function oit/tdyn (Fig.la-b) and as a func- 
tion of the position angle of the perturber 8(t) (Fig.lc-d) 
without the effects of the damped modes. We normal- 
ize to the total background energy of the primary system 
Eback = vGM'^/Rt {v = -0.78 for a King Wq = 3 model). 

For a given value of the pericenter, the peak of the re- 
sponse occurs after the closest approach. The faster the 
perturber, the larger the value of the position angle and 
the more distant the perturber location when the response 
of the primary system reaches its peak. The strength of re- 
sponse increases as the velocity of the perturber decreases. 
In Figure 2 we have plotted the maximum value of Ep as a 
function of the velocity of the perturber for the two values 
of p investigated. The response amplitude is large when 
the characteristic angular frequency, V/p, is in low-order 
resonance with primary orbital frequencies. 

Figure 3 and Figure 4 show the same plots shown in 
Figure 1 and Figure 2 but with the effects of the discrete 
weakly damped modes taken into account. The maximum 
values of Ep are about two orders of magnitude larger than 
those reached when the effects of damped modes are not 
included while the dependence of the peak of Ep on the 
velocity of the perturber and on the pericenter of its orbit 
are similar to each other in the two cases. The inclusion 
of the damped modes is critical to the persistence and the 
strength of the perturbation. 

A simple empirical function of the form K/V" can be 
used to obtain a satisfactory fit of the curves shown in Fig- 
ure 2 and Figure 4. These relations can be used to predict 
ensemble properties (see §3.4 for an example). Table 2 
summarizes the values of the parameters K and a for the 
best fits for all the cases investigated and shown in Figure 
2 and Figure 4. 

Table 2 

Best fit parameters for the scaling 

of max(0.1M/m)^_Ep/_Ebacfc with the 
relative velocity of encounters 



Parameters 


K 


a 


p/Rh — 1.0, no-damped 


0.005 


1.41 


p/Rh = 2.0, no-damped 


0.002 


1.33 


p/Rh — 1.0, with damped 


1.0 


1.96 


p/Rh = 2.0, with damped 


0.5 


1.9 



The scaling of max(0.1Af/m) Ep/Etack with the relative 
velocity of encounters, V, has been fitted using the function 

f(Y) — ^ 

^ ' (v(km/s)/200)° ' 

The scalings obtained will not extend to the very low- 
velocity regime where they would lead to a divergence of 
the energy associated to the response. In the limit of very 
slow encounters we expect the energy of the response to 

^The progressive alignment between the orientation of the wake and 
in Figure 12 and discussed in §3.3 for external encounters. 



increase more slowly and eventually to converge to a con- 
stant value (see Murali & Tremaine 1998 for an investiga- 
tion of the response of galactic halos to adiabatic pertur- 
bations). 

Figure 5 illustrates the differences between the response 
with and without the damped modes by separately show- 
ing the energy in the dipole and quadrupole components 
for V = 500 km/s and p/Rh — 1.0. As expected, the dipole 
response is significantly stronger than that of quadrupole. 
The inclusion of the damped modes leads to a much longer 
and stronger dipole response, while producing no signifi- 
cant difference in the quadrupole response. Such long-lived 
features may explain the numerous cases of galaxies with 
distorted morphologies but without any apparent interact- 
ing system in their vicinities (e.g. Richter & Sancisi 1994). 
For example, for a system with mass and radius equal to 
those of our fiducial system (M = 6 x lO^M©, R = 200 
kpc), the damping time of the dipole mode (ti = 5 x 10"'^° 
yrs) is longer the Hubble time and any perturbation ex- 
cited in a dark halo during a fly-by as well as its effect on 
an embedded disk will persist well after the encounter. 

3.2. Structure of response 

Contour plots in Figures 6 and 7 describe the density 
perturbation in the orbital plane for the dipole (upper pan- 
els) and quadrupole (lower panels) responses without and 
with the damped modes, respectively. These cases have 
V = 200 km/s and p/Rh = 1.0 and show the response 
of the primary galaxy at three different points of the per- 
turber trajectory. A lag between the orientation of the 
wake and the position angle of the perturber is evident for 
all velocities. This lag is a consequence of angular mo- 
mentum transfer between the wake and the perturber and 
tends to zero only when there is no resonance affecting 
the response of the primary system Physically, reso- 
nant momentum exchange viewed at the orbit level has 
a corresponding view as a global response: the global re- 
sponse must lag the perturber in order to apply the torque 
responsible for the exchange. 

Figures 8a and 8b show the ratio, Rp/Rh, of the radial 
position of the peak of the density of the perturbation, 
Rp, to the half-mass radius, Rh, as a function of the po- 
sition angle of the perturber without the effects of weakly 
damped modes for V = 200 km/s and p/Rh = 1.0,2.0 
(Fig.8a) and for p/Rh = 1.0 and V = 200, 500, 1000 km/s 
(Fig. 8b). The evolution of Rp/Rh is only weakly depen- 
dent on the velocity of the perturber and on the pericen- 
ter of the orbit of the perturber. The perturbation is 
efficiently transmitted to the inner parts of the primary 
system and the peak of the density of the perturbation 
is located at 0.3 < (Rp/Rh) ^ 0.4 or in terms of the core 
radius, Rc, 0.4 < (Rp/Rc) < 0.5 (for almost the entire du- 
ration of the fly-by) . The weak dependence of the response 
on perturber parameters is due to the dominant effect of 
the lowest-order point modes on the structure of the dis- 
persion matrix elements. For an analogy, consider the sim- 
ilarity of bell's ring under a variety of strikes. The strong 
oscillations of Rp/Rh in the final phases of the fiy-by are 
due to numerical difficulty in determining the density peak 
as the response weakens. 

the position angle of the perturber for very slow encounters is shown 
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Fig. 1. — Run of energy in the perturbation with time and with the position angle of the perturber for fly-bys with p/Rh = 1-0 (a) and (c) 
and p/Rh = 2.0 (b) and (d). The three curves shown in each frame correspond (from the upper to the lower curve) toV = 200 km/s, V = 500 
km/s, V = 1000 km/s. Sfcacfc is the total background energy of the primary system and t^yn is the half-mass dynamical time. Pericenter 
occurs at t = 0. 




200 400 600 800 1000 

V(km/s) 



Fig. 2. — Maximum value of the energy in the perturbation, Ep, as a function of the perturber velocity V . The upper curve corresponds 
to fly-bys with p/Rh = 1.0 and the lower curve to fly-bys with p/Rh = 2.0. 



6 



E.Vesperini, Martin D. Weinberg 




Fig. 3. — Same as Fig. la-d but including the effects of weaJdy damped modes. 




^ : : : \ : : : \ : : : \ : : : 

200 400 600 800 1000 

V(km/s) 

Fig. 4. — Same as Fig. 2 but including the effects of weaJsly damped modes. 
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Fig. 6. Density distortion of the primary galaxy at tluec pliascs of tfie fiy-by witii V = 200 km/s and p/Rh = 1-0 without the effects due 
to the damped modes. Upper panels show the dipole response; lower panels show the quadrupole response. The white dot shows the position 
of the perturber. White and black regions correspond to overdensity and underdensity respectively. 
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Fig. 7. — Same as Fig. 6 but including the effects of weakly damped modes. 




Fig. 8. — Evolution of the ratio of the position of the peak of the density distortion to the half-mass radius of the primary galaxy as a 
function of the position angle of the pcrturbcr for the fly-bys with (a) V = 200 km/s and p/Rn = 1-0 (dashed line) and p/Rh = 2.0 (solid 
line) and (b) for the flybys with p/Rh = 1-0 and V = 200 km/s (solid line), V = 500 km/s (short dashed line) V = 1000 km/s (long dashed 
line), (c) and (d) same as (a) and (b) but including the effects of weakly damped modes. 
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Fig. 10. — Energy in the perturbation as a function of the orthogonal function indices. Upper panels (a-c) correspond to the fiy-by and the 
phases shown in Fig.6; lower panels (d-f) correspond to the fiy-by and the phases shown in Fig. 7. Open circles and crosses show the energy 
in the dipole response and in the quadrupole response respectively. Each curve is normalized to the maximum value of Ep{n) . 
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Figures 9a and 9b show the evolution of the maximum 
vahio of tlic density of the perturbation Sp^ax (normal- 
ized to the central density of the equilibrium model) for 
the same cases shown in Figures 8a and 8b. The maximum 
response obtains for slower fly-bys (cf. Fig. 1) and, for a 
given value of V , Spmax increases as p decreases. Figures 
8c and 8d and 9c and 9d show the same plots as in Fig- 
ures 8a-b and Figures 9a-b for fly-bys including the weakly 
damped modes. The position of the peak of the density 
is similar to that obtained in the case without the weakly 
damped modes but there are evident differences in the time 
evolution of Rp/Rh particularly at the beginning and at 
the end of the fly-by due to the persistence of the damped 
modes. The response is significantly stronger when the ef- 
fects of the damped mode are included (compare Fig.9a-b 
with Fig.9c-d) 

Figure 10 shows the energy in the perturbation for the 
first ten basis functions and indicates the distribution of 
energy with spatial scale. Plots shown in Figure 10 are for 
the same cases and phases shown in Figure 6 and Figure 
7. Figure 10 confirms what is already qualitatively evi- 
dent in Figure 6 and Figure 7: the energy in the dipole 
response is preferentially deposited on large scales and the 
quadrupole response energy deposited on relatively smaller 
scales. The typical scale represented by each index n can 
be seen from the plots of Figure 11 which show the density 
basis functions for Z = 1. 

3.3. External encounters 

External encounters, those with trajectories everywhere 
outside the primary, are more frequent in groups and 
clusters but less dramatic. For external encounters, the 
stronger dipole response corresponds to a simple shift 
of the center of mass of the system leaving the weaker 
quadrupole as dominant aspherical contribution. Here 
we briefly describe the features of externally excited re- 
sponses. 

Table 3 summarizes the characteristics of the fly-bys we 
have investigated. An external encounter has no intrinsic 
length scale and the frequency parameter Q alone now de- 
termines the structure of the response. The Table reports 
the velocity of each fly-by assuming the pericenter to be 
equal to twice the total radius of the primary system (400 
kpc) and the maximum values of the ratio of the energy 
in the perturbation to the total energy of the equilibrium 
model. This values scales as shown for other values. 



Table 3 

External fly-bys 



n 


Vp=2Rt (km/s) 


™^^(flt/p) (™/m) Ebaak 


0.1 


23 


1.3 X 10"^ 


0.5 


115 


1.3 X 10-3 


1.0 


231 


7.0 X lO-'' 


2.5 


578 


4.7 X lO-'^ 


5.0 


1150 


2.8 X 10"^ 



Vp=2R-t is the relative velocity of the encounter for 
the given value of reported in the first column for 
for p = 2Rt= 400 kpc. 



One sees that the effects of external encounters with 
less massive perturbers are weak and they can not pro- 
duce a significant deformation in the primary system. On 



the other hand, if the perturber is more massive than the 
perturbed system even external fly-bys can be very impor- 
tant in affecting the structure of a galaxy. For example, 
in order to have max{Ep / Ef,ack) — 0.01 for an encounter 
with p = 2Rt the ratio of the mass of the perturber to 
that of the perturbed must be m/M ~ 30. These en- 
counters are relevant for spiral galaxies in clusters where 
they can undergo frequent external encounters with more 
massive giant ellipticals. A recent investigation by Rubin 
ct al. (1999) has shown that about half of a sample of 
89 spirals in the Virgo cluster exhibit signs of kinematic 
disturbances in their rotation curves; according to that 
investigation disturbed galaxies are preferentially on ra- 
dial orbits likely to cross the central regions of the clusters 
where such encounters arc likely to occur. 

The flve panels of Figure 12 show the contour plots of the 
density of the response when Ep/Eback reaches its maxi- 
mum value. As described in §3.2, the response reaches 
its maximum strength at different phases of the orbit of 
the perturber depending on the value of fi: the larger the 
value of il, the more distant the perturber is from the 
pericenter of its orbit when the response is maximum. As 
the characteristic frequency of the fly-by decreases and the 
perturbation tends to the adiabatic regime, individual res- 
onances cease to be relevant in determining the structure 
of the response and the orientation of the wake aligns with 
the position of the perturber. Finally, we note that some 
features are present at small radii and these may cause 
visible distortions in an embedded disk. 

3.4. Summary and discussion: dark halos 

We have investigated the response of a galactic dark 
halo to the perturbation induced by external and internal 
fly-bys. External fly-bys are effective when the mass of 
the perturber is larger than that of the primary system 
considered. For example, perturbations excited in spirals 
during external fly-bys with more massive giant ellipti- 
cals can play an important role in altering their structure 
and kinematics. Similarly, external encounters will pro- 
duce significant structure in dwarfs in the neighborhood 
of normal spirals. 

We have shown that the perturbation excited during 
an internal fiy-by encounter can be eflaciently transmitted 
to the inner regions of the halo where it can affect the 
structure of an embedded stellar disk. The halo distor- 
tions can be significant and persistent for typical dwarf 
interlopers. This may lead to distorted morphologies like 
those observed in warped and lopsided disks. The survey 
over different velocities and pericenter distances provides 
a quantitative measure of the expected importance of the 
process investigated in different environments. 

For an overall indicator of the fraction of objects dis- 
torted as a result of fly-by encounters, we have calculated 
the mean value of the peak relative distortion for different 
environments, 

f JL fivW^dV 
(max {Q.lM/mfEp/E,^,k) = ]f{y)V^dV 

where fiV) is the distribution of the relative velocities of 
encounters which we take equal to a Gaussian with dis- 
persion a. Figure 13 shows {va.&-x.(Q.lM / m)'^ Ep / Eback) as 
a function of a for the four scaling laws listed in Table 
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Fig. 12. — Density distortion of the primary galaxy when the perturbation energy is maximum. The three frames on the left correspond to 
fly-bys with (from the upper to the lower panel) frequencies i7 = 0.1, 0.5, 1.0; the two frames on the right are for the fly-bys with fl = 2.5, 5.0. 
The straight solid line in each panel indicates the position angle of the perturber. 
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Fig. 13. — Mean value of max(0.1M/m)^Ej,/E(,ac*; assuming a Gaussian distribution of the relative velocity of encounters versus the 
dispersion, a, of the distribution. Solid lines (upper line for p/Rh = 1, lower line for p/Rh = 2.0) are for runs including the effects damped 
modes; dashed lines are for runs without the effects of damped modes. 



2. Generally, low veloeity dispersion compact groups with 
number densities n ~ 10^ Mpc^"^ and velocity dispersions 
£7 ~ 10^ km/s are dominated by low- velocity close en- 
counters which lead to the strongest response in the per- 
turbed galaxies; cluster of galaxies which arc characterized 
by higher velocity dispersions {a ~ 10"^ km/s) and lower 
number densities (n ~ 10'^ — 10^ Mpc~^) encounters are 
expected to produce less pronounced effects. 

The effects of damped modes on the strength and the 
persistence of the features excited by the perturber are sig- 
nificant. In particular the persistence of the deformations 
due to the long-lived dipole mode may explain the ori- 
gin of the numerous lopsided disks in field galaxies which 
have no apparent companion (see e.g. Richter & Sancisi 
1994, Zaritsky & Rix 1997). We have separated the re- 
sponse due to damped modes in Figure 13 to highlight 
their importance, however, true physical response includes 
both. Damped modes may be artificially suppressed in 
n-body simulations with small particle number owing to 
fluctuation-induced orbital diffusion. 

In the next section, we will focus our attention on ellip- 
tical galaxies modeled as high-concentration King models 
for which the effects of discrete modes have not been con- 
sidered because they have much shorter damping times 
(Weinberg 1993) and thus are not as relevant as they are 
for dark halos ^. 



4. PERTURBATION OF CLUSTER ELLIPTICAL GALAXIES 



Table 4 

Orbital parameters for the fly-by 
encounters with a high-concentration 
primary system 



V (km/s) 


P (kpc) 


p/Rh 




200 


7.0 


0.5 


17.9 


500 


7.0 


0.5 


44.7 


1000 


7.0 


0.5 


89.5 


200 


14.0 


1.0 


8.9 


500 


14.0 


1.0 


22.4 


1000 


14.0 


1.0 


44.7 


200 


42.0 


3.0 


3.0 


500 


42.0 


3.0 


7.5 


1000 


42.0 


3.0 


14.9 


200 


97.8 


7.0 


1.28 


500 


97.8 


7.0 


3.2 


1000 


97.8 


7.0 


6.4 



We will model the primary system by an isotropic King 
model with Wq = 7.0 which has a concentration c = 1.53 
^. We adopt typical mass, total radius and half- mass ra- 
dius values for cluster elliptical galaxies: M = lO^^M© 
and Rt = 120 kpc, Rh = 14.0 kpc. Table 4 summarizes 
the properties of our model grid (see Table 1 for parame- 
ter definitions). The overall approach and procedure oth- 
erwise follows that described in §3 with similar physical 
interpretation. §4.1 briefly presents the response diagnos- 
tics developed in §3, highlighting the differences between 
responses in the low- and high-concentration primaries. 
§4.2 explores observational diagnostics of the fly-by in- 
duced asymmetries and proposes a generalization of one 
of the standard observational asymmetry parameters to 
test the dynamical mechanism explored here. 



*In addition, their inclusion presents challenging problems connected with the stability of the numerical evaluation of the poles and the 
residues of the dispersion matrix. 

^A core-free primary increases the cost of the computation without significant changes in the response beyond the core region 
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Fig. 14. — Evolution of the perturbation energy as a function of time (t = is the time of passage of the pcrturber at the pcriccntcr of 
its orbit) (a-c) and the position angle of the pcrturber (d-f). Panels (a) and (d) are for fly-bys with V = 200 km/s, panels (b) and (c) for 
V = 500 km/s, panels (c) and (f) for V = 1000 km/s. In each panel the curves shown correspond to different values of the pericenter of the 
orbit of the perturber: p/Rf^ = 0.5 (solid line), p/Rh = 1-0 (dotted line), p/Rh = 3.0 (dashed line), p/Rh = 7.0 (long dashed line). 




Fig. 15. — Evolution of the position of the peak density distortion as a function of the position angle for the fly-bys with V = 200 km/s 
(a), V = 500 km/s (b) and V = 1000 km/s (c). In each panel the four curves correspond to different values of the pericenter of perturber: 
p/Rh = 0.5 (solid line), p/Rh = 1-0 (dotted line), p/Rh = 3.0 (dashed line), p/Rh = 7.0 (long dashed line). 
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4.1. Structure of response 

Figure 14 shows the time evolution of the perturbation 

energy. As expected (cf. §3), the innermost fly-bys in- 
duce the strongest response and, for a given value of the 
pericenter, the response becomes weaker as the velocity 
of the perturber increases. The time evolution of Ep is 
characterized by several peaks and a non-negligible ampli- 
tude even as tends to 7r/2 and the perturber moves away 
from the primary galaxy. The persistence of these features 
well past the closest passage of the perturber may explain 
the origin of distorted morphology in galaxies without any 
close visible companion. Figure 15 shows the evolution of 
Rpl R/i. The response is manifested well inside the primary 
galaxy and it affects the very inner regions of the system 
0.05 < Rp/Rh ^ 0.15 (0.2 < Rp/Rc < 0.6) even when the 
pericenter of the perturber is as large as 3Rh or 7Rh- The 
value of Rp is determined by the dominant resonances de- 
termining the structure of the response of the primary sys- 
tem. The larger dynamic range in primary profile enables 
excitation of high-order resonances between the charac- 
teristic frequency of the motion of the perturber and the 
frequencies of the orbital motion at smaller radii in the pri- 
mary system. The rectilinear trajectory considered here 
provides continuum spectrum of frequencies that peaks at 
approximately V/p. 

We found that Rp/Rc is roughly constant for all the fly- 
bys investigated. This agrees with the results obtained by 
Weinberg (1994) in his analysis of the structure of weakly 
damped modes in isotropic King models. This implies that 
more concentrated models are more efficient in transmit- 
ting to the inner regions the perturbation excited during 
a fly- by (compare Fig. 8 with Fig. 15). 

Figure 16 shows the density distortions for a fly-by with 
V = 200 and p/Rh = 1-0 at three diff'erent phases of the 
orbit of the perturber. As in previous cases, the maximum 
response occurs well inside the primary system and the 
perturbation is efficiently transmitted to the inner regions 
of the system far inside the orbital pericenter. Figure 17 
shows plots of Ep (n) for the same case and phases shown in 
Figure 16 describing the distribution of energy at different 
scales. Similar to the results in §3, the dipole response is 
dominated by large scale features while in the quadrupole 
response there is a non negligible contribution to the total 
energy of the perturbation from features at smaller scales. 

4.2. Diagnostics of asymmetry in the primary galaxy 

We employ an established asymmetry diagnostic to com- 
pare our predictions to observations and suggest a straight- 
forward generalization which will help discriminate the re- 
sponse mechanism considered here from others. We will 
focus our attention on an asymmetry parameter A for im- 
age data, defined as 

■^rot (^5 

y)\ 

^-2 E^ra ^' 

where I{x,y) and Irot{x,y) are the brightness at a given 
point {x, y) of the original image and in the image after 
a 180 degree rotation about the position of the maximum 
density of the primary system, respectively. 

®These calculations require the evaluation of the density distortion 
tional time by approximately a factor 100 



For our calculation we use the density of the system on 
the orbital plane and thus our theoretical values are not af- 
fected by any projection effect and will provide a measure 
of the real asymmetry of the plane. Observational values, 
on the other hand, are determined using the projected den- 
sity and part of the underlying asymmetry of the orbital 
plane is likely to be washed out. In order to evaluate the 
effect of using the projected density on the orbital plane 
instead of the real plane density, we have repeated some of 
the calculations^: we have found that, for a given value of 
m/M, the values of A obtained using the projected den- 
sity are smaller than those obtained using the real plane 
density by a constant factor independent of the perturber 
orbital parameters and approximately equal to 3.8; this 
implies that, using the projected density, all the values of 
A (as well as those of the parameter S introduced later in 
this section) and the numerical parameters in the scaling 
laws of the amplitude of A with the perturber's orbital pa- 
rameters discussed in this section are not altered but they 
are produced by a perturber with a mass equal to 0.38M 
(instead of 0.1 M adopted in the rest of this section). 

The A parameter was first introduced by Abraham et al. 
(1996a) to estimate the asymmetry observed in galaxies in 
the Medium Deep Survey and in the Hubble Deep Field 
and it has been recently adopted in several investigations 
on the morphology of galax;ies in different environments 
(e.g. van den Bergh et al. 1996, Naim, Ratnatunga & 
Griffiths 1997, Marlcau & Simard 1998, Conselice & Ber- 
shady 1998, Conselice & Gallagher 1999). The values of ^ 
reported in these investigations range from for the most 
symmetric systems to 0.5 for those showing the strongest 
asymmetries. A comparison of the distribution of values 
of A for a sample of local galaxies, galaxies in the Medium 
Deep Survey and in the Hubble Deep Field was reported 
by Abraham et al. (1996a, 1996b). Their results show 
that the distribution of A for galaxies in the HDF is sig- 
nificantly skewed toward larger values of A compared to 
the distribution for galaxies in the MDS which in turn is 
skewed to larger values of A relative to the distribution of 
A for a sample of local galaxies. This trend is likely to 
be due to the higher interaction and merging rates for dis- 
tant galaxies in the HDF and the MDS. In principle, a high 
value of A could also result from an asymmetric distribu- 
tion of star forming regions, however a recent investigation 
by Conselice & Bershady (1998) shows that the colors of 
many asymmetric objects in the HDF are not blue enough 
to support this hypothesis. 

In addition, we have calculated the displacement, 5, be- 
tween the position of the peak of the density of the pri- 
mary system and the position of its center of mass. The 
estimate of S is of particular interest since many observa- 
tions have provided evidence of off-center galactic nuclei 
and of significant displacement of isophotal centers (see 
e.g. Miller & Smith 1992 for a short overview of several 
observational analyses). In particular, Naim et al. (1997) 
have adopted distances between the centers of different 
isophotal curves as one of the diagnostic parameters char- 
acterizing the degree of peculiarity of galaxies in a sample 
of moderate-redshift systems and Mendes de Oliveira & 
Hickson (1994) have investigated the morphology of early- 
type galaxies in compact groups and used the evidence of 

a three-dimensional grid which leads to an increase in the computa- 
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Fig. 17. — Perturbation energy as a function of the orthogonal function indices for three phases shown in Figure 16. In each panel, open 
circles (crosses) indicate dipole (quadrupole) response. Each curve is normalized to the maximum value of Ep{n). 
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non-concentric isophotes as a diagnostic to identify per- 
turbed galaxies. 

4.3. Evolution of the asymmetry diagnostic parameters 

Figure 18 shows the time evolution of the asymmetry 
parameter A defined above (eq. 0) . Although several well- 
documented processes contribute to the observed distorted 
morphologies, the majority of the values in the distribu- 
tion of A obtained for galaxies in the MDS and in the 
HDF by Abraham et al. (1996) (see their Fig. 2) fall well 
in the range of values of A we have obtained in our anal- 
ysis. Mergers and strong interactions with more massive 
systems are likely to be responsible for the production of 
the extreme values of A (A > 0.3) observed in some galax- 
ies. We point out that the observational data refer to all 
the morphological types including spiral galaxies to which 
our results do not directly apply. 

Figure 19 shows the time evolution of S. 

Table 5 

Best fit parameters for the scaling of 

max(0.1A//m)A and max(0.1Af/m)5 (pc) 
with the relative velocity of encounters 



id. 


Ka 


aA 


Ks 


as 


p/Rh = 0.5 


0.24 


0.6 


2.5 X 10^ 


0.6 


p/Rh = 1.0 


0.19 


0.5 


1.9 X 10^ 


0.5 


p/Rh - 3.0 


0.11 


0.4 


0.8 X 10^ 


0.4 


p/Rh = 7.0 


0.03 


0.1 


0.1 X 10^ 


0.0 



the relative velocity of encounters, V, has been fitted using 
the function f{V) — 



(y(km/s)/200)° 



As expected, the trend with V and p is similar to that ob- 
tained for the asymmetry parameter A, with slow and close 
encounters being those capable of producing the largest 
values of S. Even fly-bys with large values of p are capa- 
ble of affecting the inner regions of the primary galaxies 
and lead to a significant displacement of the position of 
the density peak from the center of mass of the system. 

The dependence of both asymmetry parameters on V 
and p is summarized in Figure 20 which shows the maxi- 
mum values taken by A and S" as a function of p/Rh for 
the three values of V considered. We have fit the scaling 
of max0.1(M/m)^ and max0.1(M/m)S' with the relative 
velocity of the encounter for different values of p with a 
power-law K/ {V{ km s^^)/200)" and we report the values 
obtained for the cases investigated in Table 5. These fits 
are not valid in the limit of adiabatic encounters where the 
strength of the perturbation and of the diagnostic param- 
eters tend to a constant value as described previously (see 
§3.1). 

4.4. Radial dependence of the asymmetry parameter A 

Fly-by interactions described here produce structure on 
well defined scales (cf. Fig. 16), however, the sum in 
equation (^ washes out the signature of our predicted 
response by averaging over regions with little distortion. 
To improve sensitivity, we propose an incomplete form of 
the asymmetry parameter A{r) which restricts the sum 
in equation (|^) to pixels enclosed by a circle of radius r 
from the position of the peak of density Rp. Figure 21 
shows the plot of A{t,r) (calculated for m/M = 0.1) for 



the case with V ~ 500 km/s and p/Rh — 0.5. The appear- 
ance of A{t, r) is qualitatively similar for the other cases 
investigated. The profile of A at any given time has a 
maximum at 0.3 < r/R^ ^ 1.0 with a subsequent decrease 
which is steeper near the absolute maximum of A{t, r) and 
shallower far from this maximum. The value of A always 
tends to a constant for large values of r. In Figure 22a we 
compare the radial profile of A for the three fiy-bys with 
p/Rh = 1.0 and V = 200, 500, 1000 km/s. The three 
curves show the run of A{r) at peak and each curve is nor- 
malized to the maximum value of A(r), Amax- The radial 
location of Amax and the shape of A{r) do not depend 
significantly on the velocity of the perturber in these half- 
mass units. Figure 22b compares A{r)/Amax for all the 
fly-bys with V = 500 km/s and with the different values of 
p/Rh investigated. The location of the peak of A{r)/Amax 
depends very slightly on p/Rh while there is a marked dif- 
ference between the slope of A{r)/Amax after the peak for 
the fiy-bys with small pericenters {p/Rh < 1.0) and those 
with larger values of p/Rh {p/Rh > 3.0). A distinct peak 
in the radial profile of A is present only at some phases of 
the encounter and this could be difficult to capture in ob- 
servations. On the other hand the decrease at small values 
of r is a feature common to all the phases of the orbit of 
the perturber and it should be an observable signature. 



4.5. Summary and discussion 

In this section, we have studied the response of a high- 
concentration spherical stellar system to the perturbation 
induced by a fly-by. We represent the outer profiles of 
cluster ellipticals by a high-concentration King model; we 
choose a King profile rather than a core-free model for 
numerical convenience. A deformation in the structure of 
the primary galaxy is directly observable and we explored 
quantitative measures of the asymmetry produced during 
the fly-by encounter. In the end, we focussed our attention 
on the asymmetry parameter A, introduced by Abraham 
et al. (1996a, 1996b) and measured for a number of local 
galaxies and in galaxies in the MDS and in the HDF fields, 
and on the shift S, between the position of the peak of the 
density of the system and its center of mass. We present 
the time evolution, the scaling with the perturber orbit 
and the radial dependence of the asymmetry to facilitate 
estimates for a wide variety of environments. 

The values of A caused by fiy-bys range from A w to 
A « 0.2 and this same range of values of A is found in 
observations of numerous local galaxies and distant galax- 
ies in the MDS and in the HDF fields (Abraham et al. 
1996a, 1996b). Extreme values of A observed for some 
galaxies in the HDF field {A > 0.2) are probably caused 
by mergers. The quantities A and S reach their maximum 
values when the perturber is close to the periccntcr of its 
orbit. As evident in Figure 18, the distortions may be 
long-lived; in some cases the perturber can be quite dis- 
tant from the primary system and therefore we can expect 
some non-negligible values of A and S even without any 
close companion in the vicinity of the primary. 

Our results show that the values of A and S will be 
determined by both the distribution of relative encounter 
velocities and the density of the environment. In short, 
richer groups and clusters with lower velocity dispersions 
will lead to high values of A and S. Figure 18 shows that 
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Fig. 18. — Evolution of the asymmetry parameter A (see eq.[7|) as a function of the position angle for the fly-bys with V = 200 km/s 
(a), V = 500 km/s (b) and V = 1000 km/s (c). In each panel the four curves correspond to different values of the pericenter of perturber: 
p/Rh = 0-5 (solid line), p/Rh = 1-0 (dotted line), p/Ru = 3.0 (dashed line), p/Rh = 7.0 (long dashed line). 
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Fig. 19. — Evolution of the distance, S, between the position of the density peak and its center of mass as a function of the position angle 
of the perturber for the fly-bys with V = 200 km/s (a), V = 500 km/s (b) and V = 1000 km/s (c). In each panel the four curves correspond 
to different values of the pericenter of perturber: p/Rh = 0.5 (solid line), p/Rh = 1.0 (dotted line), p/Rh = 3.0 (dashed line), p/Rh = 7.0 
(long dashed line). 
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Fig. 21. — Plot of the asymmetry parameter, A, as a function of time and radius of the area (centered at the position of the peak of density) 
considered for its calculation for the fly-by with V = 500 km/s and p/Rh = 0.5. 
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Fig. 22. — (a) Radial profile of A (normalized to its maximum value) for fly-bys with p/Rh = 1-0 and V = 200 km/s (solid line), V = 500 
km/s (dotted line), V = 1000 km/s (dashed line); each curve correspond to the radial profile of A when the total A is maximum, (b) Same 
as (a) but for fiy-by with V = 500 km/s and p/Rh = 0.5 (solid line), p/Rh = 1-0 (dotted line), p/Rh = 3.0 (dashed line), p/Rh = 7.0 (long 
dashed line). 



disturbances persist until O > 1.4 or roughly t ^ h x 10^ 
years after the pericenter passage. In order to provide a 
quantitative indication of the dependence of the response 
on the environmental conditions, we have calculated the 
probability that a galaxy has an encounter in the time 
interval r such that max A > 0.1 assuming that the per- 
turber has a mass m = O.IM (or m = 0.38M if the pro- 
jected density is used for the calculation of A). For our 
fiducial system t = 5 x 10^ yrs and m = IQ^^AdQ, the 
probability P of having such an encounter is given by 



mTp^{V)VT 



(8) 



where n is the galaxy number density of the environment 
and piy) is the pericenter distance for which max A = 0.1 
during an encounter with velocity V. The function p{V) 
has been determined using the parametrization reported 
in Table 5 for the values of V available from our numerical 
investigation and then fitting these points with a function 
of the form K/Vl^. We obtain 



P{n,V) = 



10-3n(Mpc"^) 
(y(km s-^)/200)i-44 



(9) 



and describe this function in Figure 23. This quantifies 
our finding that galaxies in environments characterized 
by high density and low-velocity dispersion, like compact 
groups, will be more likely to show signs of past or recent 
interactions in agreement with the results of a number of 
investigations (see e.g. Zepf & Whitmore 1993, Mendes 
de Oliveira & Hickson 1994). These observational studies 
argue that the fraction of distorted morphologies in com- 
pact groups is larger than the fraction observed in clus- 
ters (which have similar number densities but higher ve- 
locity dispersions) and in sample of field galaxies (which 
are likely to suffer low-velocity encounters but in a low- 
density environment). Assuming that young galaxies were 



located in dense environments and characterized by lower, 
previrialized velocity dispersions, our conclusions are con- 
sistent with the results of Abraham et al. (1996b) who 
have shown that the distribution of A for galaxies in the 
HDF is skewed toward high values compared to that for 
galaxies in the MDS which, in turn, are more asymmetric 
than local field galaxies. 

An additional useful diagnostic is the ratio of the per- 
turber to the primary masses necessary to produce a given 
value of max A. This is easily obtained from the fits dis- 
cussed in §4.3 and is given by 



' TO N 



0.1 



max j4 



y(km/s) 
200 



(10) 



where the values of Ka and a a depend on p/Rh (cf. Table 
5). The four frames in Figure 24 show the contour plots of 
(m/A'/)maxA in the plane max A— (as we discussed at the 
beginning of this section, if the projected density is used 
for the calculation of A the values of to/M are larger than 
those reported in Fig. 24 by a factor 3.8). Any environ- 
ment producing low-velocity, close encounters can lead to 
significant asymmetries (max A > 0.1) in the primary sys- 
tem even by low-mass perturbers (to/M ~ 0.02 — 0.1). For 
example, consider a system with a mass and a radius equal 
to those of our fiducial system (M = 10^^ ji = ^20 
kpc). An encounter with a perturber with m/M = 0.05 
with relative velocity V — 200 km/s leads to max A = 0.1 
ii p/Rh = 1.0 {p ~ 14 kpc for our fiducial system). Low- 
mass interlopers with pericenters close to the primary cen- 
ter could lead to significant asymmetries in cluster ellipti- 
cals as well. 

Fly-by encounters yield the following signature in our 
radially-dependent asymmetry parameter, A{r): 1) a de- 
crease at small radii, 2) a peak around 0.3 < r/Rh ^ 1.0 
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Fig. 23. — Loci of constant probability P{n, V) for an encounter in the past t = 10t<jj,„ such that max A > 0.1. The vertical and horizontal 
axes describe the number density of the environment and relative velocity of encounters. 
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p/Rh=1.0 
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Fig. 24. — Loci of constant relative pertuiber mass necessary to have a maximum value of the asymmetry parameter equal to max A during 
an encounter with relative velocity V. Pericenter distances label each frame. 
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and 3) a convergence to a constant value at large radii. 
Although we have not predicted the signature of A{r) for 
scenarios other than fly-bys, we believe an observational 
investigation would provide a test of the mechanism pro- 
posed here and would be a general probe of asymmetric 
structure. The computation of A(r) is relatively straight- 
forward with the data already available. 

5. CONCLUSIONS 

We have explored the structure and the persistence of 
the features in dark halos and cluster ellipticals excited by 
an unbound (or fly-by) interaction over a range of differ- 
ent galactic environments. The main conclusions of our 
analysis are as follows: 

1. Fly- by encounters which penetrate the outer regions 
of a galaxy can lead to asymmetric features well in- 
side the half-mass radius. The distortion of a dark 
halo may affect its disk (e.g. Weinberg 1998b). 

2. In dark halos the lowest-order damped modes yield 
strong features which persist well after the per- 
turber's passage and could account for distorted 
morphologies in galaxies without obvious close com- 
panions. Such damped modes are a natural part of 
the coUisionless response but are difficult to resolve 
in N-body simulations with modest particle number 
{N < 10^) due to noise. 

3. The strength of the response depends both on the 
velocity of the perturber and on the pericenter of 
its orbit; low- velocity, close encounters produce the 
strongest perturbations. 

4. Completely external fly-by encounters will produce 
a strong response if the perturber is more massive 
than the primary. For example such encounters can 
produce significant effects on spiral galaxies in clus- 
ters where they are likely to interact with more mas- 



ioi'inal galaxies. 

5. A proposed generalization of the asymmetry param- 
eter A (Abraham et al. 1996a), provides a diagnostic 
of the radial structure of the asymmetry and a test 
of the mechanism explored here in particular. Our 
predicted range of A is similar to that observed for 
the majority of galaxies in the Medium Deep Sur- 
vey, in the Hubble Deep Field and for local galaxies. 
We quantify the range of high-density, low-velocity 
dispersion environments expected to produce signif- 
icant distortions by fly-by encounters. Indeed, com- 
pact groups have been shown to have a significant 
fraction of distorted galaxies. 

Additional discussion is provided in §3.4 for dark halos and 
in §4.5 for cluster ellipticals. 

Future work will explore the properties of the response 
for cuspy galaxies; the magnitude of the resonant coupling 
to the inner galaxy in core-free profiles would nicely com- 
plement the present work. We plan to further investigate 
the effects of the halo distortions on an embedded stel- 
lar and gaseous disk, by extending our current study to 
explore the response of the disk and to quantify the in- 
duced deformations. This will permit direct comparisons 



with the observed distortions in the morphology of disks 
as well as to study possible effects of these distortions on 
their kinematics. In particular, the effects of even small de- 
formations in the morphologies of disks could significantly 
affect their kinematics possibly giving rise to the observed 
scatter in the TuUy-Fisher relation as suggested by Franx 
& de Zeeuw (1992). 

We thank Eric Linder for comments on the manuscript. 
This work described here was supported in part by NSF 
AST-9529328. 
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APPENDIX 

DETAILS OF THE DERIVATION OF THE MATRIX EQUATION 

Our results are based on an explicit solution for the response of a spherical stellar system to the perturbation induced 
by a perturber with mass m on rectilinear trajectory with pericenter p. Specifically, we solve the coupled linearized 
coUisionless Boltzmann equation for the distribution function of the primary stellar system and Poisson's equation (cf. 
§2) 

dt aw 91 ai 9w ' ^ ' 

V2$i = 47rG/9i, (2) 

where the subscript denotes the equilibrium quantities and the subscript 1 denotes the first order perturbation of the 
distribution function /, and Hamiltonian H. The perturbed potential, $i, is the sum of the tidal potential due to the 
perturber, <i>p, and the gravitational potential of the response to the perturbation, '^resp- After performing the Laplace 
transform in time, the Fourier transform in the angle variables and solving for /i, equation (|l|) becomes 

(Ai) 

where we denote the Laplace transform of a quantity y{t) by y[s) and the Fourier transform of a quantity y(w) by y-y 
with 7 = (71, 72, 73) indicating the vector of integer indices for the discrete Fourier expansion in the angle variables. The 
frequencies associated to the angle variables are denoted by w = dH^ / dl. 

Following the approach adopted in Weinberg (1989), we will use a potential-density biorthogonal basis, u\, d\ to expand 
the perturbed quantities $1 and pi. The pairs u', d\ are solutions of the Poisson equation. The basis we have adopted 
is obtained by numerical solution of the associated Sturm-Liouville problem in which the lowest order basis functions are 
tailored to the background equilibrium model following Weinberg (1999). The solution to the coupled equation, which 
determines the response of the halo to the perturbation, then reduces to the solution of an algebraic integral equation for 
the coefficients of the expansion. The subsections b elow sketch the main steps in the derivation of the algebraic equation 
for the coefficients of the expansion from equation ( |AlD (see e.g. Weinberg 1989 for additional details). 

Expansion of the perturbation potential in action-angle variables 

We begin by expanding the perturbing potential, $p, in terms of spherical harmonics Yim{0, 4>) (where 9 and denotes 
the angular coordinates inside the primary stellar system). Up to quadrupole order, the potential of a point mass perturber 
can be written as 



$p = Gm 




(A2) 



where &{t) = arctan(Ut/p) is the position angle of the perturber defined as the angle between the axis from the center of 
the system to the pericentric position of the perturber and the axis from the center to the position at time t. We define 
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the usual quantities r< = min(r, i?), r> = max(r, i?) where r and R are the radial coordinate inside the stellar system 
and the distance between the center of the system and the perturber, respectively. 

Following Tremaine & Weinberg (1984), the radial part of the perturbing potential is expanded in terms of the basis 
function u\ and each resulting term is expressed in the form u\{r)Yijn{d, (j)) a Fourier series in the angle variables w. 
We finally obtain: 



7=-oo 



- -V27273^730 + \ ^27273*^732 + \ l^e '"t^27273^' 



/27273O73-2 I ^27273 



^— ^ (-e *®Vi7273(573i + e*®Vi7273(573_i) X 



71 

17273 



where 



and 



17273 27r / 73 \ ' 



(A3) 



(A4) 



(A5) 

The angle /3 is the inclination of the orbital plane to the equatorial plane and ri^^^^ {(3) are the components of the rotation 
matrices (see e.g. Edmonds 1960). 

Derivation of matrix equation 
Similarly, we can express the perturbing and the response potential in the form 

* = E ^7273 (/?) E < ii)^'^;^/^'^- (A6) 

1,1 J 

For the perturbing potential, each coefficient a^^^ {t) can be expressed as the product of a time-independent coefficient, 

b\^^, and a function of time depending only on I and 73, gi^^{t). The functional form of gi-y^it) can be easily derived for 

any value of / and 73 using equation ([A^). For example, for r < R, I — 2 and 73 — 2, one finds 322 — cos^ 0(i)e~^'®^*^. 
Equation (|A1[) now takes the form 



/71 



s + i^ ■ uj 



E ^'7273 E "'73 (^)^r7273 e'^-" + E ^'7273 iP) E ^/733'73 (^)^77273 e'^-^ 



(A7) 



- P ~ resp 

Pi ^ Pi 



For self-consistency, the response density, which is obtained by integrating equation (A7) over the velocity coordinates, 
presp^ must be equal to the density pf determined by the Poisson's equation with $1 = (j)^'^''^. Imposing the condition 

(A8) 
)wing 

(A9) 
(AlO) 

(All) 
(A12) 



multiplying both sides of equation( A8) by Yi*^^u\ and integrating over the spatial coordinates, we obtain the following 
equation for the coefficients 

«j73(^)=E<3(^)[«/73(^)+^/735(^)]- 

j 

Equivalently, this may be written as a matrix equation: 

a,73(s) = D,:^^(s)M,73(s)[b,735(s)] 
where the elements of the matrix M/73 (s) are defined as 

^^3(^)- E / dIdw-^-^;^l,,_,3l^7,73X-4^^3X7-^ 



s + i-f ■ U) 



and 

Equations 
(no externa 



Di73 = I - M 



'73- 



or (AlO) completely determines the response to the perturbation induced by the perturber. For b = 
perturbation), one obtains the following nonlinear eigenvalue problem 



^'73 



(A13) 



The eigenvalues of this problem arc the frequencies of the point modes of the stellar system: frequencies with Re{s) > 
correspond to unstable growing modes while for stable damped modes Re{s) < 0. 
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Inverse Laplace transform 

Solution of equation ( [A9| ) describes the evolution in frequency space. In order to get the time evolution of the pertur- 
bation we need to perform the inverse Laplace transform of equation (A9). This yields the following expression for the 
coefficients Sii^ysit) 



it) 



E 



'73' 



(A14) 



where, to simplify the notation, we have defined 



kj _ 



dl 



'7273^;-72-73 /7273 ' 



(A15) 



Equation (A14) does not take into account the zeros of det D(s) which are the point modes of the primary system. Since 
we will consider only systems known to be stable, there are no modes with Re{s) > 0, while modes with Re{s) < will 
be damped and die away with time. Nevertheless, as shown in Weinberg (1994), for King models some modes can be very 
weakly damped and persist long after their excitation. When the effects of damped modes are included. Equation (A14) 
is valid for t — > cx) while for finite time it has an additional term for each such mode. Altogether we find 



*i73 



71 :72 
ft 



J <J —OO J <J 



gi^,{t )e-'^' dt 



bi 



If 3 ■ 



kj 



(A16) 



where Sd denote the frequencies of t he da mped modes. 

Numerical evaluation of equation (A16) requires truncating the infinite sum over 71 (the sums over the indices 72 and 
73 range from —I to I; see Tremaine & Weinberg 1984). After checking the convergence of the solution, the final value 
adopted is "fimax = 6 for the low-concentration King model studied in §3 and "fimax — 10 for the high-concentration King 
model considered in §4. Similarly, the sum over the radial basis must be truncated. Since the the adopted biorthogonal 
functions are tailored to the equilibrium model Umax = 10 is sufficient for the convergence. All the integrals have been 
calculated by Romberg's method with the number of points chosen so to guarantee a maximum error of « 10^^. 



